Dissipative quantum systems: from two to many atoms 
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We study the dynamics of bosonic atoms in a double well potential under the influence of dissipation. The 
main effect of dissipation is to destroy quantum coherence and to drive the system towards a unique steady 
state. We study how the atom-atom interaction affects the decoherence process. We use a systematic approach 
considering different atomic densities. We show that, for two atoms, the interaction already strongly suppresses 
decoherence: a phenomenon we refer to as "interaction impeded decoherence". For many atoms, thanks to the 
increased complexity of the system, the nature of the decoherence process is dramatically altered giving rise to 
an algebraic instead of exponential decay. 
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I. INTRODUCTION 



Every device, may it be classical or quantum, operates sur- 
rounded by an environment which affects to a certain degree 
its properties. Understanding the interplay between the differ- 
ent physical processes at work and the dissipative processes 
induced by the environment is thus of paramount importance. 
This is of particular interest for the development of technolo- 
gies based on quantum mechanics as the environment can, for 
example, destroy the quantum coherence between different 
parts of the device and thus completely spoil its functional- 
ities. Various quantum systems affected by dissipative pro- 
cesses have been studied in solid-state physics fllHjl, atomic 
and molecular physics [5] and quantum optics (a [J. How- 
ever, at the moment, a complete understanding of the physics 
of strongly correlated quantum many-body systems coupled 
to an environment is lacking. In recent years, due to rapid 
technological progress, cold atom experimentalists have been 
able to build quantum emulators, systems mimicking with a 
high degree of precision modeled Hamiltonians. For this rea- 
son, these cold atoms have become an ideal test-bed to study 
quantum dynamical systems. In fact, the theoretical study of 
ultracold atoms paved the way to the first understanding of the 
interplay between interactions and dissipation. Some notable 
works are (i) the discovery of Zeno-like dynamics due to local 
losses of atoms in a double well geometry ifsl- TToll or in optical 
lattices (ii) the existence of augmented stability against 
dissipation for weakly II 1211 and strongly 11131 1 1 4fl interacting 
systems and (iii) the emergence of a dynamical phase tran- 
sition between the condensed and thermal steady states lfl5ll . 
Also of great importance are works detailing methods to engi- 
neer dissipative processes generating dark states with precise 
quantum properties lHa - COB . 

In this work we will study a quantum well in a systematic 
way by investigating different boson densities. The system 
under study is a collection of ultracold atoms in a double well 
potential under the influence of a dissipative process. As a 
warm-up, we first consider a double well loaded with only 
two atoms. This exercise will help us gain a deeper analyt- 



ical understanding of the combined effect of interaction and 
dissipation, and show that this interplay leads to a marked 
slowing down of the decoherence process. We then study a 
system with many atoms and demonstrate that the increased 
complexity of the system results in the emergence of a regime 
where coherence decays algebraically. Hence, for large atom 
numbers and interaction strength the system is particularly re- 
sistant to dissipation. 

This paper is organized as follows. In Section [II] we intro- 
duce the model and the various relevant parameters. In Sec- 
tion|III]we first study in detail the case of two atoms in a dou- 
ble well before investigating in|IV]the dynamics for the same 
double well loaded this time with a large number of atoms. In 
Section[V]we draw our conclusions. 



II. THE MODEL 

The properties of an ensemble of bosonic atoms in a double 
well can be described, in the single band approximation, by 
the Hamiltonian 



H 



(« 



J b[b 2 + 



ftj(hj - 1). (1) 



Here bj (St) is a bosonic operator annihilating (creating) an 



J 3 ^"3< 

atom at site j and hj = SjSj counts the number of atoms 
at site j. Dealing with a double well, j can only take the 
values 1, 2. The parameter J is the hopping amplitude, and 
U the interaction strength. As the ratio U/J increases, the 
phase coherence between the two sites, C — (b\b2 + SJjSi), 
decreases. To represent a state vector, we use the Fock basis 
\l,m) where / and m label the number of atoms on sites 1 and 
2, respectively. 

As this ultracold gas of bosons is under the influence of a 
dissipative process, the system can be described by a density 
matrix p which evolves following the master equation 



df' 



H,p +iD{p) 



(2) 



2 



where D(p) is the dissipator modeled by 



(3) 



3=1,2 



with coupling strength 7. It has been shown in fl 1 3L 12111 that 
this dissipator, written in Lindblad form 1I22I l23ll . can be ex- 
perimentally realized, for example, in the context of ultracold 
atoms exchanging energy with a red-detuned optical lattice. 
Let us note that the considered dissipator leaves the Fock 
states invariant. Here we study the situation where the sys- 
tem is first prepared in the ground state of the Hamiltonian. 
At time t = the dissipation is switched on and the evolution 
of the system under the influence of both the Hamiltonian and 
the dissipator is investigated. 



in. TWO ATOMS 

Having established the model and main notations in the pre- 
vious section we now attempt to gain a deeper understanding 
of the physics of this system. The best starting point is to 
study a double well loaded with two atoms as this limit can be 
solved analytically. We show here that the system converges 
towards a unique steady state and that the relaxation process 
slows down as the interaction strength increases. 

As a reminder, we first provide here a brief summary of the 
physics of the Hamiltonian system without dissipation. The 
dimension of the Hilbert space is three and we choose as a ba- 
sis the states |2, 0), |1, 1), |0, 2) in this order. The Hamiltonian 
in this basis is thus given by 



H = 



U -V2J 
-V2J ~V2J 
-V2J U 



(4) 



The ground state of this Hamiltonian, which we choose 
as the initial state for the time evolution, has an energy 

Eg = 



U — \JU 2 + 16J 2 ) and is given by tpc 



a(l, (U + VU 2 + 16 J 2 ) /y/SJ, 1), where a is the normal- 
ization constant a = 1/^2+ (U + y/U 2 + 16J 2 ) 2 /8J 2 . 
Note that for large U / J the ground state wavefunction tends 
towards (V2J/U, 1 - 2J 2 /U 2 , V2J/U) with energy E G w 

u ■ 

With these results in mind, we can now address the fol- 
lowing question: how does the system evolve under the ef- 
fect of the dissipation considered? Using ©, we compute the 
time evolution of the nine elements of the density matrix p. 
The evolution of the six complex independent elements 112411 
is given by 





= iV2J (p 1)2 


- P2,l) 




(5) 




= iV2J (p 2 ,i 


~ P\,2 + 


Pa.i - Pi,o) 




hpo,o 


= iV2J (pifl 


~ P0,l) 






hp2,i 


= iv2J(pi,i 


- P2,2 - 


P2,o) - iUp2,i 






= iV%J (pi.o 


~ P2,l) - 






hpi,o 


= iV2J (p 2 ,o 


+ Po,o ~ 


Pis) + iUpifi 





Given that we prepare the system in the ground state of the 
Hamiltonian, the initial state is symmetric under the exchange 
of the first and third element. For this reason we only focus 
on the symmetric subspace of the density matrix evolution. 
Eq. (0 is a linear equation and its symmetric eigenvalues A are 
the roots of the polynomial p(X) = 48(hj) 2 J 2 + [4(?i7) 3 + 
6«7 J 2 + Ah~fU 2 ] A+ [9(ft 7 ) 2 + 16 J 2 + U 2 ]X 2 + 6hjX 3 + A 4 , 
and the steady state, with eigenvalue A5 = 0, corresponds to 
the density matrix 



PS 



1 
1 
1 



(6) 



which is the completely mixed state (or the highest entropy 
state). The roots of p(X) are in general complex numbers with 

Re(X) < 0. 



A. Dynamical transitions 

Studying the roots of p(X), we identify dynamical transi- 
tions which arise when tuning U/J and Tvy/J transforms a 
real into a complex eigenvalue (or the converse). Dynami- 
cally this implies that a simple exponential decay, the hall- 
mark of real eigenvalues, can be transformed into an exponen- 
tial decay accompanied by oscillations (due to the imaginary 
part). The oscillations have their origin in the Hamiltonian 
evolution. The presence of dynamical transitions becomes 
particularly interesting if one changes the value of the differ- 
ent parameters, for example U or 7, in time, allowing for the 
modulation of the dynamics from an overdamped to an under- 
damped regime. The various dynamical scenarios are depicted 
in Fig. Q]for different values of tvy/J. For tvy/J < 6.3, as 
show in Fig.[TJa,b), two eigenvalues are real and two are com- 
plex conjugates. For 6.3 < tvy/ J < 6.8, Fig.[Ttc,d), there are 
three regimes as a function of U/J: in the first regime, there 
are two complex and two real A, in the second regime all the A 
are real, and in the last regime again two eigenvalues are real 
and two are complex. For larger values of the dissipative cou- 
pling, fir// J > 6.8, there are two regimes: for small Uj J, the 
decays are only exponential, while for large U/J two eigen- 
values are real and two are complex (see Fig.[TIe,f)). 

For any value of tvy/ J, if U /J is large enough the imag- 
inary part of two of the eigenvalues will be proportional to 
±/7/ J indicating the strong influence of the Hamiltonian part 
of the dynamics. 

The change in dynamics is exemplified in Fig. |2] where 
tvy / J = 6.5. In this figure, we plot the absolute value of the 
coherence versus time in a log-log form which highlight the 
change of sign of the coherence for U / J = 0.1 (blue line). 
We observe that, as expected, for both U/J — 0.01 (blue 
solid line) and U/J = 20 (red dot-dashed line) the evolu- 
tion of the coherence is oscillatory while for U/J — 1 (green 
dashed line) there are no oscillations but only a sum of expo- 
nential decays. The inset highlights the oscillatory behavior 
for U/J = 20. 
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FIG. 1: Real (a,c,e) and imaginary (b,d,f) part of the eigenvalues 
A (divided by J) versus the interaction U/J for hy/J — 4 (a,b), 
hy/J = 6.5 (c,d) and h-y/J = 7 (e,f). 
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FIG. 2: Absolute value of the coherence versus time in a log-log plot 
for hy/J = 6.5 and for interaction U/J — 0.1 (blue solid line), 
U/J=l (green dashed line), U/ J — 20 (red dot-dashed line). The 
inset shows the coherence versus time in a linear plot to highlight the 
oscillating behavior for U/J = 20. 



B. Quantum Zeno effect versus interaction impeded 
decoherence 

The long time dynamics of the system are governed by its 
slowest decaying state, those whose corresponding eigenval- 
ues have the smallest (in modulus) real part. For this small 
system with only five symmetric states, the relevant ones are 
the steady state ps and the next slowest decaying state pe 
with corresponding A^. While the steady state is indepen- 
dent of U/J, pe and its decay rate can change considerably 
with U/J. To gain a deeper understanding of how the interac- 
tion affects the dynamics, we first present the non-interacting 
case U — and then cover the large interaction limit where 
U/J > 1 and U/h-y > 1. 



1. Quantum Xeno effect 

For the non-interacting case, U — 0, the equation for 
the eigenvalues p(X) takes the simpler form p(X)u=o = 
(ft 7 +A) {48fr 7 J 2 + [4(/i 7 ) 2 + 16J 2 ]A + 5ft 7 A 2 + A 3 }. The 
slowest decaying eigenvalue at large h-y/J is, in this regime, 



-12 



try 



(7) 



This is a clear manifestation of the Quantum Zeno effect: in- 
creasing the dissipative coupling 7 suppresses the decay rate 
as the continuous density measurement prevents the system 
from evolving. An experimental observation of this effect was 
reported in II25II . 



2. Interaction impeded decoherence 

When the interaction is the dominant energy scale (for large 
U / J and large U /tvy), we analytically find that 



A, 



A/ 



-12 



ft 7 J 2 
U 2 



(8) 



Similarly to the Quantum Zeno regime, here too the decay is 
strongly suppressed. However, as the physical mechanism be- 
hind this slowing down is different from the Quantum Zeno 
case, the dependence of A/ on the system parameters is dis- 
tinct from Xqz- Remarkably, in this regime, the presence of 
interactions renders the system more robust against dissipa- 
tion, as it increases the size of the energy gap that needs to 
be overcome to populate higher states. As reaching the steady 
ps requires a large amount of energy, of the order of U, the 
dissipative process, heating up the system, will need a longer 
time to provide this energy at large U. Consequently, we refer 
to this phenomenon as "interaction impeded decoherence". 

The cross-over between the two regimes is clearly illus- 
trated in Fig. [3j a) where, for small interactions, A decreases 
as 7 increases, while, for large interactions, A and 7 increase 
together as indicated by the green arrows. 

For higher clarity we have also shown in Fig. |3jb,c) 
the value of, respectively, (Xe/J) x {Ivy /J) an d 
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FIG. 3: (a) Slowest decaying even eigenvalue Xe/J versus inter- 
action U/J (blue diamonds) for different fry/ J = 10,20,30,40. 
The blue diamonds indicate the exact numerical values while the red- 
dashed line represents the approximate solution from ( |13t . (b) Xe/J 
times fry/ J versus interaction U/ J. (c) Xe/J divided by fry/J ver- 
sus interaction U/J. The arrows indicate increasing 7 values. 



{\e / J) / {hf / J) . The collapse of the curve in Fig. |3b) 
confirms the emergence of a Quantum Zeno behavior while 
the collapse of the curve in Fig. [3c) signals the emergence of 
the "interaction impeded coherence" behavior. 



C. Large dissipation or large interaction limit 

We gain deeper insight into the large interaction U ^S> J 
and/or large dissipation regime fry ^> J via an adiabatic elim- 
ination technique II 14 12611 . The main idea behind this method 
is that at long times the dynamics is dominated by the diag- 
onal elements of p which dictate the full quantum dynamics. 
In fact, their evolution relates to the hopping parameter J, the 
smallest parameter in the system, and connects the diagonal 
and the off-diagonal elements. However, the off-diagonal el- 
ements evolve very quickly, with a dephasing due to the in- 
teraction U and/or a fast exponential decay due to 7. Hence, 
the evolution of the off-diagonal terms can be integrated out 
and we are left with a coarse-grained evolution for the whole 
system. 

More precisely let us take, for example, the evolution of 
po,i from (0. This is given by 



Po,i(i) 



-(i¥+7)t 



Po,i(0) 



j^ e -(^+7)t 

h 



is much smaller and decays even faster. Also as U and 7 are 
large, we can integrate (O by parts to obtain lUil I2H 



po,i(t) ft v2J 



Pi,i - Po,o 
U — ifry 



Likewise we have 

Pi,a(t) ft 

P0,2 ft 



fK T U + ifi-y , 



V2J — ^ 7 xn (p 2 ,2 



U 2 + (fry) 2 



Po,o) 



Pl.l) 



(10) 



0. 



(11) 



We now inject the information from ([Tot and (flU into the time 
evolution of the diagonal elements and, as anticipated earlier, 
the diagonal values of p are determined by three independent 
equations. This results in 



-77 G°0,o + Pi,i + P2,2) = (12a) 
at 

■4 (Po,o ~ P2.2) ft ^ (po,o ~ p2,2) (12b) 
at n 

d . . Xe / n \ 

{P0,0 - 2/51,1 + P2,2) ft ~T \P0,0 - 2/01,1 + P2,2) ■ 



dt 



(12c) 



Eq. d!2a| ) shows consistently that ps is the steady state, dl2b| ) 
refers to the slow decaying anti-symmetric solution, which de- 
cays atarate Xo/h — — 4J 2 7/([/ 2 +(fry) 2 ) butis notrelevant 
for our dynamics, and the last equation, (1 12c) , represents the 
symmetric slowest decaying state with a rate 



h ~ C/ 2 + (fry) 2 ' 



(13) 



From ( fT3] l, we derive both the Quantum Zeno behavior of Q 
and the "interaction impeded decoherence" behavior (|8) in 
their respective regimes (small and large interactions limits). 

Using ([Tol l and ( fTTT i. we compute, at any given time, the co- 
herence of the system knowing only the terms on the diagonal. 
In particular, the coherence of the slowest decaying state A e 
is given by 



C E oc 



JU 



U 2 + (fry) 2 



(14) 



which is proportional to U when it is small compared to fry, 
and inversely proportional to U when it is the dominant energy 
scale. Thus the large interaction strengths affect the coherence 
in two ways: U lowers the coherence of the slowest decaying 
state, but at the same time reduces the decay rate pE- 



(9) 



IV. MANY ATOMS 



Po,o — Po,2) dr. We now wa nt to understand how the dynamics is altered by 
the presence of many atoms. As for the two atom case the 
Here it should be noted that po,2, being far from the diagonal, interaction was found to play a crucial role in the dynamics, 



5 



we can expect that the presence of many interacting atoms will 
affect the dynamics even more dramatically. We first present 
a perturbative study at short times highlighting the subtle role 
played by interactions, we then cover in detail the strongly 
interacting regime where coherence decays algebraically with 
time. 



strength, the short time evolution is slowed down as U op- 
poses the effect of 7. This constitutes a first glimpse into the 
role played by interactions. In the following, we derive an ex- 
pression for the decay of the coherence C at large times and 
for a large number of atoms. This expression provides a much 
deeper understanding of the emerging physics. 



•'Interaction impeded decoherence": first seeds 



B. Eigenvalues problem 



To illustrate how interaction influences decoherence, we 
study perturbatively the time evolution of the coherence. We 
first show analytically that for a non-interacting system the 
coherence decays following a simple exponential. To prove 
this statement it is sufficient to show that the evolution of the 
operator C = b\l>2 + b\b\ follows 



dt 



C 



H,C 



1 

n 

= -ic. 



■ 7 ©(C) 



(15) 



This expression is derived from the fact that T>(C) = C and 



that, for C/ = 0, 



H,C 



0. In the non-interacting limit, the 



evolution of the coherence is then simply described by 

C(t) u=0 = e-^C(0). 



(16) 



For non-zero interaction strength, and in particular for large 
values of U, the decoherence follows a different law even for 
short times. To demonstrate this statement, we analyze the 
evolution of the density matrix p up to third order in dt 



p(t + dt) = p(t) + 'p(t)dt + f^-dt 2 



6 



(17) 



where each superimposed dot stands for a time derivative. As 
we use the ground state of the Hamiltonian as the initial state, 

we have that H, pa =0 and hence can write 



V(V(p G ))-- H,V(p G ) 



PG = V{p G ) 
PG 

- l - \h,V(D(pg)) 



(18) 



H,V(M 
1 r 



K 2 



H, 



H,V{p G ) 



The coherence is only dependent on the real part of the off- 
diagonal terms of the time-evolved density matrix, hence, for 
U I J 1, we derive 



Cull) = [ l-jdt+-^-dt 2 



7 3 - l{U/hf 



dt 3 j C(0). 

(19) 

If U = 0, this expression would be the Taylor expansion of 
an exponential decay. However, due to the finite interaction 



A very good first insight is provided by the study of the 
eigenvalues A of the evolution of the density matrix p. For the 
general case with N atoms, the size of the density matrix p is 
(N + l) 2 . The elements of the density matrix pi :Tn are such 
that p = J2i m Pl,m\h N — I) (m, N — m| where N is the total 
number of atoms. Rewriting the matrix p in a vector form p, 
we compute at each time t the vector p. This vector is given 
by the sum over all p a , each associated with a complex decay 
constant X a , 



Pit) 



e Kt/h 



(20) 



here c a is the weight of each p a at time t — 0. The real 
parts of the different A Q are plotted in Fig. [4] Two aspects of 
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FIG. 4: Real part of the eigenvalues Re (X a ) as a function of the 
number of atoms N for fey/ J = 1 and U/J = 20. The red-dashed 
line represent an algebraic fit that scales as 1/N 2 . 

this figure are worth noting: (i) the slowest decaying state is 
gapped from all others (this gap depends on the number of par- 
ticles as 1/N 2 ) and (ii) for large N an almost continuous band 
of slowly decaying states arises. The first point implies that 
past a certain critical time the decay will be fully exponential 
and solely due to the slowest decaying state ps (note that the 
corresponding A e decreases as 1 /N 2 as shown by the linear 



6 



fit in Fig.|4]i. This state, as expected, is mainly populated on 
the diagonal, see Fig. [5] and these diagonal elements follow 
a particular distribution. We observe that the matrix elements 




FIG. 5: Typical plot of the density matrix p for the slowest decaying 
symmetric state. This particular case is for TV = 30, U/J = W and 

hj/J = l. 



corresponding to large density imbalance between the left and 
right sites, for example po,o or Pn,n, are much more pop- 
ulated than those with balanced configurations, for example 
Pn/2,n/2- For a large interaction strength, states with strong 
imbalance are highly energetic and thus populating them re- 
quires the dissipation to provide a tremendous amount of en- 
ergy. Moreover, in order to populate the state p nn , with n 
large, it is necessary to first populate state p n -\. n -\ which in 
turns necessitates a long time. Thus, the different configura- 
tions can be arranged into an energetical hierarchy. 

The second aspect, the presence of a band of eigenvalues, 
paves the way to the emergence of a new decay regime. This 
interesting behavior arises as similar exponential decays are 
summed over. The emergence of this regime is shown in Fig. [6] 
where coherence decay for various values of U / J and TV is 
plotted. On the left panels (a,c,e), we show how the slower 
decaying states decay in time, while on the right panels (b,d,f) 
are plotted the weighted coherence of each eigenstate given by 
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FIG. 6: Left column: decay of the coherence C(t) versus time 
for different values of TV and U/J. Right column: decay of the 
weighted coherence of each eigenstate C a (t) versus time. The in- 
set of Fig.[6fe) presents a log-log of the coherence and highlight the 
power-law regime. The parameters are TV = 20 and U / J = 2 for 
(a) and (b), TV = 5 and U/J = 20 for (c) and (d), and TV = 20 and 
U/J = 20 for (e) and (f), while fir// J = 1 for all the subplots. 



C. Emergence of an algebraic relaxation regime 



To provide a better understanding of this emerging dy- 
namical behavior and to pinpoint the algebraic nature of this 
regime, we first proceed with a coarse-graining of the evolu- 
tion using adiabatic elimination and then the resulting equa- 
tions are further approximated within a continuum limit. The 
end product of this procedure is a continuum classical diffu- 
sion equation with a non-uniform diffusion function. The so- 
lution of this diffusion equation is then used to demonstrate 
the existence of the algebraic regime [ 14]. 



C a (t) 



tr 



(CPc 



(21) 



For small U/J or TV the exponential decay of the slowest 
decaying state pE takes place earlier as the interaction is 
not strong enough to impede decoherence (small interaction 
U / J = 2 in Fig.|6ja,b)) or there are not enough states to form 
a band (small atom number TV = 5 in Fig.|6jc,d)). When both 
the interaction and atom number are sufficiently large, for ex- 
ample U/J — 20 and TV = 20 as in Fig. [6je,f), a power-law 
decay regime emerges over a large time. This regime is clearly 
shown on the log-log plot in the inset of Fig.|6je). 



1. Adiabatic elimination 

As for the two-atom problem described in section III, we 
study the long time full quantum dynamics by only consid- 
ering the diagonal elements of p. The adiabatic elimination, 
in the regime of large dissipation Jvy/ J and interaction UJ J, 
thus reduces the size of the system to (TV + 1). The system of 
equations describing the dynamics is derived by generalizing 



7 



([Tol l to the limit [//Try > 1 to obtain 



Pn,n+l(t) 



Jy/(n + l)(N-n) 



U 2 {N -2n-l) 2 
x (U(N - 2n - 1) - ihj) 

x (/?n+l,n+l ~ Pn,n)- (22) 

This generalization is also applied to the master equation (0 
and gives 

N 

9tPo,q = ^ N _ x y Oi,i - Po.o) (23a) 
_ (n + l)(N-n) 

O r Pn,n — -777 — \Pn+l,n+l ~ Pn,n) \^->") 



9 T pN 



(7Y-2n-l) 2 
n(iV-n + l) 
+ (iV- 2n+ l) 2 
N 



[Pn— X,n— 1 Pn,n) 



N 



{pN-i,N-i - Pn,n) (23c) 



(N- l) 2 

for 71 ^ {0, TV}, JV even and r = 2(J 2 /U 2 )jt. These last 
equations are used to accurately compute the dynamics and 
give access to the coherence of the system 

^ ^ 2J(n+l)(N-n) f 

° ~ / J TTl nT — 7^ — \Pn+l,n+l — Pn,n)- 



U(N -271-1) 



(24) 



Continuum limit 



By performing the large limit, the discrete master equa- 
tion ( 1231 is mapped onto a classical diffusion equation. The 
configuration space, indexed by 71 is mapped to the coordi- 
nate x = n/N — 1/2 € [—1/2, 1/2] which is a continuous 
variable in the limit A*" —> 00. The boundaries of this system 
then become x = ±1/2 corresponding to the strongly imbal- 
anced configurations of the double well, n = and n = N, 
whereas the center, x = 0, corresponds to the balanced config- 
uration n = N/2. The diagonal elements of the density ma- 
trix p are connected to a continuum probability density distri- 
butions p(x, t) via the relation N p n ^ n {r) = p(x, r). Within 
this limit, the normalization of the probability density distri- 

1 /2 

bution is given by J_y 2 p(x)dx = 1 which is equivalent to 
tr p = 1. The initial state corresponds to p(x, r = 0) = 5{x), 
and the steady state is then given by the uniform distribution 
p(x, t = 00) = 1 which represents the totally mixed state p$. 
To derive the differential equation for the probability density 
function p(x, r), we rewrite (123b) as 

(x + dx)(l — x) 



d T p{x, t) 



(1 - 2x - dx) 2 
x(l — x + dx) 
(1 - 2x + dxf 



{p{x + dx, t) — p(x, t)) 



(p(x- dx,r) -p{x,r)) 
(25) 

where dx = 1/N. Expanding the right-hand side of d25l l, to 
second order in dx, gives the diffusion equation 



d TN p(x,t n ) 



d x [D(x)d x p(x, tn)] 



(26) 



where tn = t/N 2 . The diffusion function D{x) = — 1 
varies substantially as a function of x: for the balanced config- 
uration x — 0, D(x) diverges while for strongly imbalanced 
configurations x (=> ±1/2 it vanishes. Thus, if the initial state 
is near the origin, which is the case when the evolution begins 
from the ground state, we observe a very fast initial diffu- 
sion; when the probability density distribution approaches the 
boundaries the diffusion drastically slows down. 

For short times, such that p(x,tm) is not affected by the 
boundaries, the diffusion equation ( |26] i is solved by using a 
scaling ansatz 



1 



p(x,r N ) = — /(£) 



(27) 



with £ = x/t 1 ^. A scaling solution exists for v — 1/4 and 
tat <C 1 giving the differential equation for / 



if + (c 4 - 2)/' + e 3 / = 0. 

This expression is solved analytically giving 

/(0 ex exp(-£ 4 /4). 



(28) 



(29) 



The diffusion for the probability density function p(x, tn) is 
thus, for short rescaled time tn <C 1, given by 



P(X,T N ) 



V2 



1 



r(i/4) r v 4 



exp (-x 4 /4-nv) 



(30) 



where the constant r(() is the gamma function with argu- 
ment (. Eq. (f30b describes an anomalous diffusion process 
as (x 2 ) oc ^/tn which, since this formulation is valid for rjv 
small, and as for normal diffusion (a; 2 ) oc tn, ( f30l > is an ex- 
ample of super-diffusion. Hence, the divergence in D{x) at 
x « leads to a highly accelerated initial diffusion. 

Regarding the coherence C, it is possible to obtain an ana- 
lytical expression from the continuum limit of (l24l : 



C_ 

N 



J 

UN 
J 



1 ' 2 x 2 



1/4 



-1/2 

r(3/4) 



d x p(x, r N )dx 

r(3/4) 



1 



UN2T{\/A)^ ^(1/4)^' 



(31) 



(32) 



This expression is in excellent agreement with the numerical 
simulations (see Fig. |7). Hence, we demonstrated here how 
the power-law regime emerges for a large number of atoms 
and strong interactions. Fig. [7] also shows clearly that if the 
interaction is not strong enough (black continuous curve for 
U/J= 0.1) the power-law regime does not emerge. It also 
shows that the algebraic region will increase when more atoms 
are in the system. In fact the power-law region is larger for 
large N, for example N = 60 (blue squares), and it is smaller 
for fewer atoms, for example N — 8 (light blue stars). From 
Fig. |2l we also notice that for smaller 7 the algebraic region 
begins at later times and is thus reduced (see green diamonds 
curve for N = 50, U/J = 10 and h/J = 0.1). The fact 
that the time-evolution of the coherence is independent of the 
interaction strength is particular to the coherence in this prob- 
lem. 
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FIG. 7: (color online) Coherence versus time for N = 60, U / J — 20 
and hy/J = 1 (blue squares), for N = 50, U/J = 10 and fvy/J = 
0.1 (green diamonds), for N = 8, C//J = 20 and fry/ J = 5 (light 
blue stars) and for TV = 40, U/J = 0.1 and h'y/J = 1 (black 
continuous line). The red dashed line represents the analytical result 

of CD- 

V. CONCLUSIONS 

To summarize we have shown here how the interaction be- 
tween atoms can affect decoherence in a quantum system. We 



analyzed a double well setup by first considering only a sys- 
tem of two atoms and by subsequently extending our study 
to a system of many atoms. For two atoms, we discovered 
that various dynamical behaviors emerge and depend on the 
ratios between U/J and U/fvy. We also showed that strong 
interactions suppress strongly the detrimental effects due to 
dissipation. We refer to this process as "interaction impeded 
decoherence". When the double well is loaded with many 
atoms, we found that the decoherence process is slowed down 
by interactions and that the nature of the decay is altered and 
follows an algebraic law. 

The experimental realization of this model is possible us- 
ing cold atoms trapped in optical lattices. One of the main 
challenges would be to use a system that can be well de- 
scribed by a single band model despite the large interaction 
strength. This could be achieved using a very light species 
like Lithium whose interaction can be tuned using a Feshbach 
resonance 11271 12811 . An alternate route would be to use two 
different kinds of bosonic atoms trapped in a single well po- 
tential and to control the effective tunneling and interaction 
respectively via two-photon combined microwave and radio- 
frequency pulses, and a Feshbach resonance as in ll29ll . 

We thank P. Barmettler, H.P. Breuer, J. Dalibard, J.-P. Eck- 
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sions. We acknowledge ANR (FAMOUS), SNSF (Division II, 
MaNEP), CIFAR, NSERC of Canada and the DARPA-OLE 
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